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Introduction 

We intend to comment on some of those aspects of the theory of differential 
equations which we think are clarified (for us, at least) by means of the syn- 
thetic method. By this, we understand that the objects under consideration 
are seen as objects in one sufficiently rich category (model for SDG), allowing 
us, for instance, to work with nilpotent numbers, say d E R with (P = 0; 
but the setting should also permit the formation of function spaces, so that 
some of the methods of functional analysis, become available, in particular, 
the theory of distributions. 

The specific topics we treat are generalities on vector fields and the solu- 
tions of corresponding first- and second-order ordinary differential equations; 
and also some partial differential equations, which can be seen in this light, 
the wave- and heat-equation on some simple spaces, like the line R. For these 
equations, distribution theory is not just a tool, but is rather the essence of 
the matter, since what develops through time, is a distribution (of heat, say), 
which, as stressed by Lawvere, is an extensive quantity, and as such behaves 
covariantly, unlike density functions (which behave contravariantly); and the 
distributions may have no density function, in particular in the setting of 
model for SDG where all functions are smooth. 

When we consider these partial differential equations, we shall follow an 
old practice and sometimes denote derivative d/dt with respect to "time" by 
a dot, /, whereas differential operators with respect to space variables are 
denoted df/dx, /', A(/), etc. 

We want to thank Henrik Stetk^r for useful conversations on the topic 
of distributions. 



1 Generalities on actions 



Recall that an action of a set (object) D on a set (object) M is a map 
X : D X M ^ M, and a homomorphism of actions (M, X) — > (X, Y) is a 
map f : M ^ N with /(X(rf, m)) = F(rf, /(m)) for all m G M and d e D. 

The category of actions by a set D form a topos; we shall in particular 
be interested in the exponent formation in this topos, when the action in the 
exponent is invertible. An action X : D x M ^ M is called invertible, if 
for each d & D, X{d,—) : M — >■ M is invertible. In this case, the exponent 
{N, yY^'^^ may be described as N'^ equipped with the following action by 
D: an element d E D acts on /5 : M — > X by "conjugation": 

l3^Ydol3o{Xa)~\ 

where Y^ denotes Y{d, —): N N, and similarly for X^- 

In the applications below, D is the usual set of square zero elements in R. 
It is a pointed object, pointed by G Z), and the actions X : D x M ^ M we 
consider, are pointed actions in the sense that X{0,m) = m for all m G M, 
or equivalently, Xq : M —>■ M is the identity map on M. A pointed action, 
in this situation, is the same thing as a vector field on M, cf. ||12|| . 

In the above situation, if X and Y are pointed actions, then so is the 
exponent described. The pointed actions likewise form a topos, and the 
exponent described is then also the exponent in the category of pointed 
actions; cf. 0. 

For the case of vector fields seen as actions by i^, we want to describe 
the "streamlines" generated by a vector field in abstract action-theoretic 
terms; this is going to involve a certain "universal" action (i?. A): R is an 
"infinitesimally open subset" of R, i.e., whenever x & R then x + d E R for 
every d E D. The main examples of such subsets are R itself, the non-negative 
numbers -R>o, open intervals, and the set D^o of all nilpotent elements of 
the number line. The action A is the vector field d/dx, meaning the map 
D X R ^ R given by {d, t) i— * d + t. (So it is not to be confused with 
the Laplace operatot A, to be considered later.) The main property to 
be assumed is that the individual A^'s are homomorphisms of D-actions 
(which is a commutativity requirement); the structure of R could probably 
be derived from this, but we shall be content with assuming that R is an 
additively written monoid, and that D (1 R (with the of D also being the 
zero of the monoid). 
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First, if (M, X) is a set with an action, a homomorphism / : {R, A) — 
(M, X) is to be thought of as a particular solution of the differential equation 
given by X, with initial value /(O), or as a "streamline" for the vector field 
X, starting in /(O). One wants, however, also to include dependence on 
initial value into the notion of solution, and so one is led to consider maps 

F : Rx M ^ M, 

satisfying at least F{d,m) = X{d,m) for aX\ d E D and m G M; we shall 
consider and compare the following further conditions (universally quantified 
over all d G -D, t, s G -R, m e M): 

F(A(c?,t),m))=X(rf,F(t,m)); (1) 

this is the main one, the two following conditions are included for systematic 
reasons only: 

F(A(rf,t),m)) = F(t,X(rf,m)), (2) 

F(t,X(rf,m)) = X(rf,F(t,m)) (3) 
Finally, one may consider the following equation 

F(t + s,m) = F(t,F(s,m)). (4) 

Writing X^ for the map X{d, — ) : M — > M, and similarly for F, condition 
([^) may be rewritten as 

-FA(<i,t) = XdO Ft 

The others may be rewritten in a similar way. For instance (0) may be 
rewritten as 

Ft+s = FtoF, 

Equation (|ip expresses that, for each fixed m G M, the map F{—,m) : R — »• 
M is a homomorphism (and thus, by virtue of F{0,m) = m, a "solution 
with initial value m"). Writing the action of D in terms of the symbol ■, we 
may write it F{d ■ t,m) = d ■ F{t,m). Equation (Q) expresses a certain bi- 
homogeneity condition of F, F{d ■t,m) = F{t, d-m) ; (§) says that for fixed 
t & R, F{t, — ) : M — > M is an endomorphism of D-actions, F{t, d ■ m) = 
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d ■ F{t, m). Finally (Q) is the usual condition for action af a monoid on a set 
M. Clearly, it implies all the others. 

Let X be a vector field on M, thought of as a first-order differential 
equation. We say that the map F : R x M ^ M is a. complete solution or 
simply a solution if = Xd and F satisfies (|lD. A solution in this sense does 
not satisfy the other conditions (^-(H), but it does, provided that M satisfies 
a certain axiom (refiecting, synthetically, validity of the uniqueness assertion 
for solutions of differential equations on M). — The axiom in question is the 
following 

Uniqueness property for M: 

If X is a D- action on M , and f,g:R^M are homomorphisms of actions, 
with f{0) = g{0), then f = g. 

Note that the validity of the axiom, for a given M, depends on the choice 
of R, A. For instance, we shall prove below that it holds for any microlinear 
M if i? is taken to be Doo (and A = d/dx). 

Proposition 1 Let X he a vector field on M and assume that M satisfies 
the uniqueness axiom. Then any solution F : Rx M ^ M of the differential 
equation X satisfies properties and Furthermore, if R is a monoid 
(under +) then F also satisfies 

Proof. Since the proofs are quite similar, we shall do only (^. Fix 
m E M and s G i? and define the couple of functions /, g : R —>■ M hj the 
formulas 

fit) = F,+,(m) 
g{t) = Ft o Fs{m) 

We have to check that / and g are homomorphisms of D-actions, i.e., they 
satisfy (|ip. Let us do this for the first 

f{t + d) = F(^t+d)+sim) 
= Fd+(t+s){m) 
= FdO Ft+s{m) 
= Xdofit). 

The proof that g is a, homomorphism is similar. Thus, the equality of the 
two expressions follows from the uniqueness property assumed for M. 
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Recall that a vector field X on M is called integrable if there exists a 
solution F : RxM ^ M. If we assume the uniqueness property, the equation 
(|^) holds; if further the commutative monoid structure + on i? actually is a 
group structure, then (^) implies that the action is invertible, with X^a as 
X^^ (in fact F^d = F^^)- Of course, both the uniqueness property and the 
question whether or not the vector field X is integrable, depends on which 
R is considered. In particular, we shall say that X is formally integrable or 
has a formal solution if X is integrable for R = Doo (which is a group under 
addition). For the case of M = i?", this amounts to integration by formal 
power series, whence the terminology. 

Theorem 2 The uniqueness property holds for any microlinear object, (for 
R = Doo). Furthermore, every vector field on a microlinear object is formally 
integrable. Thus, every vector field on a microlinear object has a unique 
formal solution. 

Proof. We need to recall some infinitesimal objects from the literature on 
SDG, cf. e.g. [O. Besides D <Z R, consisting of d E R with ci^ = 0, we 



have C the ra-fold product of D with itself. It has the subobject 
D{n) C D"' consisting of those n-tuples {di, . . . , dn) where di ■ dj = for all 
i,j. There is also the object Dn C R consisting of 6 & R with 5"+^ = 0; D^o 
is the union of all the D^s. If {di, . . . , dn) G D"-, then rfi + . . . + (i„ G Dn- 

— Now, let M be a microlinear object, and X a vector field on it. We 
first recall that if di,d2 G D have the property that di + d2 G D, then 
Xd^ o Xd2 = Xrf^+rfj. (For microlinear objects perceive D{2) to be a pushout 
over {0} of the two inclusions D -0(2), and clearly both expressions given 
agree if either di = or d2 = 0.) In particular, Xd^ and commute. But 
more generally. 

Lemma 3 If X is a vector field on a microlinear object and di,d2 G D, the 
maps Xd^ and Xd^ commute. 



Proof. This is a consequence of the theory of Lie brackets, cf. e.g. 3.2.2, 
namely [X, X] = 0. 

Likewise 

Lemma 4 If X is a vector field on a microlinear object and di, . . . ,dn E D 
are such that di + . . . + dn = 0, then 
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Xd^ o . . . o Xd„ = 1m 

(= the identity map on M). In particular, {Xd)~^ = X_d. 

Proof. We first prove that R, and hence any microhnear object, perceives 
Dn to be the orbit space of under the action of the symmetric group Sn in 
n letters: Assume that p : —>■ R coequahzes the action, i.e. is symmetric 
in the n arguments. By the basic axiom of SDG, p may be written in the 
form 

p{di, ...,dn)= ^Q^^ 

QC{l,...,n} 

for unique oq's in R (where d'^ denotes Hieg'^i)- We claim that aq = a^(Q) 
for every vr G S„. Indeed, 

J2aQdQ = niY.aQd'^) 
Q Q 

since p is symmetric. But 

Q Q 

By comparing coefficients and using uniqueness of coefficients, we conclude 
= a7r(Q), and this shows that p is (the restriction to of) a symmet- 
ric polynomial i?" R. By Newton's theorem (which holds internally), p 
is a polynomial in the elementary symmetric polynomials cij. Recall that 
ai{di, . . . , dn) = di + . . . + dn'. and each cxj, when restricted to D", is a 
function of Ui, since df = 0; e.g. 

0-2(^1,(^2) = '^didj = ^{di + . . . + dnf = ^((Xi(rfi, . . .,dn)f. 

Now consider, for fixed m G M, the map p : D"' — > M given by 
[di, . . . , dn) I— > Xdj^ o . . . o Xd„{m). By Lemma |^, this map is invariant under 
the symmetric group S„ (recall that this group is generated by transposi- 
tions), so there is a unique : Dn M such that 

(f){di + . . . + dn) = Xd^ o . . . o Xd„ (m). 

So if rfi + . . . + d„ = 0, Xrf, o . . . o Xd„{m) = 0(0) = 0(0 + . . . + 0) = 
Xq o . . . o Xo(m) = m. This proves the Lemma. 
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We can now prove the Theorem. We need to define Ft : M M when 
t G -Doo- Assume for instance that t G -D„. By microhnearity of M, M 
perceives Dn to be the orbit space of under the action of S„ (see the 
proof of Lemma ID , via the map (rfi, . . . , dn) ^ di + ... + dn, so we are forced 
to define Ft = X^^ o . . . Xd„ if F is to extend X and to satisfy @. The fact 
that this is well defined independently of the choice of n and the choice of 
di, . . . ,dn that add up to t follows from Lemma ^ 

As a particular case of special importance, we consider a linear vector 
field on a microlinear and Euclidean i?-module V. To say that the vector 
field is linear is to say that its principal-part formation \^ ^ ^ is a linear 
map, A, say. We have then the following version of a classical result: 

Proposition 5 Let a linear vector field on a microlinear Euclidean R-module 
V he given by the linear map A : V ^ V . Then the unique formal solution 
of the corresponding differential equation, i.e., the equation F{t) = A{F{t)) 
with initial position v, is the map D^o x V V given by 

{t,v)^e'-\v), (5) 

where the right hand side here means the sum of the following "series" (which 
has only finitely many non-vanishing terms, since t is assumed nilpotent): 

v + tAiv) + '^-A\v) + '^^A%v) + ... 



Here of course A^(f ) means A{A{v)), etc. 

Proof. We have to prove that F{t) = A{F{t)). We calculate the left 
hand side by differentiating the series term by term (there are only finitely 
many non-zero terms): 

Hv) + ^_-A\v) + -"-^A^v) + ... = A{v + t. A{v) + '^_-A'{v) + ...) 

using linearity of A. But this is just A applied to F{t). 

There is an analogous result for second order differential equations of the 
form F {t) = A{F(t)) (with A linear); the proof is similar and we omit it: 
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Proposition 6 The formal solution of this second order differential equation 
F= AF, with initial position v and initial speed w, is given by 

Fit) =v + t-w+ ^-Aiv) + ^A(«;) + ^^A\v) + ^A^H + .... 



2 Exponent vector fields 

In this section, we sliow tliat solutions of an exponent vector field may be 
obtained by conjugating solutions of the vector fields that make up the ex- 
ponent. Furthermore, this method of conjugation is equivalent (under some 
conditions) to the method of change of variables, widely used to solve differ- 
ential equations. 

Theorem 7 Assume that (M, X) and {N,Y) are vector fields having solu- 
tions F : R X M —>■ M and G : R x N ^ N, respectively, and assume that 
all Ft are invertible. Then a solution H : R x M ^ M of the exponent 
{N,Y)^'^'^^ is obtained as the map 

H : Rx N^^ ^ 
given by conjugation: Ht[(3) = GfO /3 o Ff^. 

Proof. This is purely formal. For f3 G A^^^, we have 

{Y^UH.m = Y,oH,{(3)oXj' 

= Y,oG,opoFi'oX^' 
= G,+,o(3oF^_,\ 

= H,+m, 

where in the third step we used the equation (|lD for G and F, in the form 

Gd+t = YdoGt, respectively Fd+t = XdoFt, 

together with invertibility of F^ for all s and invertibility of Xd- 

A similar argument gives that if each of (@)-(§) holds for both F and G, 
then the corresponding property holds for H. 



8 



In most applications, the invertibility of the Ft will be secured by sub- 
traction on R, with Ff^ — F_t. 

Recall that an J?-module V is called Euclidean if the canonical map a : 
V xV ^ given by a{u, v){d) = u + d-v is invertible; the composite of 
with projection to the second factor, —>■ V x V —>■ V is called principal 
part formation, li X : V ^ is a vector field on a Euclidean module V, 
we may compose it with principal part formation to get a (not necessarily 
linear) map ^ : V ^ V , called the principal part of the vector field X] it is 
thus characterized by the formula 

X{v){d)^v + d-i{v). 

Recall also that if /3 : M — > is any map into a Euclidean i?-module, and 
X is a vector field on M, then the directional derivative Dx{P) of P along X 
is the composite 

where the last map is principal part formation. Using function theoretic 
notation, Dx{P) is characterized by validity of the equation 

/3(X(m, d)) = P{m) + d ■ Dx{P){m), 

for all d e -D, m e M. 

When M itself is a Euclidean R module, and X has principal part ^, we 
usually write D^{j3) instead of Dx{l3). 

Proposition 8 Assume that Xi, X2 are vector fields on Mi, M2, respec- 
tively, and that H : Mi M2 is a homomorphism (i.e., it preserves the 
D-action). Let V he a Euclidean R-module. Then for any u : M2 — > V , 

DxAuoH)^Dx,{u)oH. 

Proof. This is a straightforward computation: 

u{X2{H{m), d)) = u{H{m)) + d ■ Dx,{H{m)); 
on the other hand 

u{X2{H{m), d)) = u{H{Xi{m, d))) = u{H{m)) + d-Dx,{uo H){m). 
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By comparing these two expressions we obtain the conclusion of the Propo- 
sition. 



For any object A^, let us consider its "zero vector field" Z , i.e., is the 
identity map on A^, for all d. For a vector field X on an object M, we then 
also have the "vertical" vector field Z x X on x M. 

If we have a complete solution F:i?xM— s>Mofa vector field X on M, 
we may consider the map F : Rx M ^ Rx M given by (t, m) ^— (t, F{t, m)) 

Proposition 9 The map F thus described is an automorphism of the vector 
field Z X X onRx M. 

Proof. By a straightforward diagram chase, one sees that this is a restate- 
ment of d^). 

We now consider solutions F : R x V ^ V for such vector fields, so 
equation (|I]) holds: X^ oFt = Ft+d- In terms of principal parts, this equation 
may be rewritten as 



Using directional derivatives, we can give a more familiar expression to 
the vector field (lODE) Y-^ considered above on the object X*^, when the 
base iV is a microlinear Euclidean i?- module V , and the exponent M is 
mocrolinear. In fact, letting rj be the principal part of the vector field Y on 
N = V,we have, for u e me M, de D (recall that {Xd)'^ = X_d) 



(at the third equality sign, a cancellation of d - d took place in the last term) 
In other words, the principal part of Y^ is 6 : M ^ V given by 




Ft = D^{Ft). 



(6) 



{Y^Uu){m) 



Ydouo X_d{m) 

M((X_d(m)) -t- d ■ r/(M(X_d(m))) 
u{m) — d ■ Dx{u){m) + d ■ rj{u{m)) 
u{m) + d ■ [—Dx{u){m) + rj{u{m))\ 



6{m) = ri{u{m)) — D x{u){m) . 
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Recalling that the lODE corresponding to a vector field X on a Euclidean 
i?-module V may be written as x = C,{x) where ^ is the principal part of X. 
In these terms, the above equation may be rewritten (leaving out the m, and 
modulo some obvious abuse of notation) as 

u = 7]{u) - Dx{u), 

or still, recalling that (— ) is "derivative with respect to time 
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— + Dx{u) = r]{u). 

This is a PDE of first order "in time" . 

The following may be seen as a generalization of (P), and is a form of the 
chain rule. We consider a vector field X on M, with solution F : RxM —>■ M. 
Let U : R X M ^ V he any function with values in a Euclidean i?-module. 

Proposition 10 Under these circumstances, we have 

^f/(t,F,(m)) = ^(t,F,(m)) + {Dz^xU){t, Ft{m)) 

for all t e R, m e M. 

Proof. Since F is a solution of X, Ft+d = XdoFt, and so for any t,t' E R 
{Z X X)d{t' , Ft{m)) = (t', Ft_|_(i(m). Therefore, by definition of directional 
derivative, 

f/(t',Fi+rf(m)) = U{t',Ft{m))+d-{Dz>cxU){t',F,{m)). 
Putting t' = t + d, we thus have 
Uit + d, Ft+dim)) = U{t + d, Ft{m)) + d ■ {Dzy,xU){t + d, Ft{m)) 

= Uit + d, Ft{m)) + d ■ {Dz^xU){t, Ft{m)) 

by a standard cancellation of two rf's, after Taylor expansion. Expanding the 
first term, we may continue: 

= f/(t, Ft{m)) + d-^{t, Ft{m)) + d ■ {Dz^xU){t, Ft{m)). 
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On the other hand, 

U{t + d, Ft+d{m)) = U{t, Ft{m)) + d ■ ^U{t, Ft{m)); 

comparing these two expressions gives the result. 

The method of change of variables has been used extensively to solve 
differential equations. We shall prove that our method for solving the ex- 
ponential differential equation , where X is an integrable vector field on 
M, Y an integrable vector field on a Euclidean i?-module, and where R is 
symmetric with respect to the origin (if t & R, then —t G R), may be seen 
as an application of the method of change of variables. We let r] : V ^ V 
denote the principal part of Y, as before. Let F : R x M ^ M he the 
assumed solution of X, and let F : R x M ^ R x M he the map 

F{t,m) = {t,F{-t,m)) 

Then F (which represents the change of variables r = t, ji — F{—t,m)) is 
invertible. 

Theorem 11 ("Change of variables"). If u : R x M —>■ V is a particular 
solution ofY^, or, equivalently, of 

^ + D^{u)^r,{u), (7) 

then the unique map U : R x M ^ V given as the composite 

RxM ^^X' RxM 



is a particular solution ofY^, or, equivalently, of 



dU 

-^ = n(u), (8) 



and vice versa. 



Proof. Since u{t,m) — U{t, F_t{m)), we have 
Ou dU 

— {t,m) = -U{t,F_t{m)) = —{t,F_t{m)) - Dz^xU{t,F_t{m)), 
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by the chain rule, Proposition [1^. On the other hand, F is an automorphism 
of the vector field Z x X, by Proposition ||, and so, by construction of F and 
Proposition ^, 

DzMu) = DzMU o F) = {Dz^xU) o F. 

Therefore, 

du 

= + Dzxxiu) - g{u) 

dU 

= — (t, fi) - DzMU)it, fi) + DzMU){t, li) - v{U{t, /i)), 
where /i = F_t{m), i.e., U is solution of 

proving the theorem (the vice versa part follows because F is invertible). 

Example. Let D be the set of elements of square zero in R, as usual. 
It carries a vector field, namely the map e : D x D ^ D given by {d, 6) 
{1 + d) ■ 6. It is easy to see that this vector field is integrable, with complete 
solution E : R X D ^ D given by {t, 5) i— >■ e* ■ 6. Now consider the tangent 
vector bundle on M. The zero vector field Z on M is certainly integrable, 
and so we have by the theorem a complete integral for the vector field Z'^ 
on the tangent bundle. We describe the integral explicitly (this then also 
describes the vector field, by restriction): it is the map R x — > 
given by {t,(3) i— > [c? i— >■ /3(e~* ■ d)]. — The vector field on obtained this 
way is, except for the sign, the Liouville vector field, cf. p[, IX.2. 



3 Generalities on distributions 

We want to apply parts of the general theory of ordinary differential equations 
to some of the basic equations of mathematical physics, the wave- and heat- 
equations. This takes us by necessity to the realm of distributions. Not 
primarily as a technique, but because of the nature of these equations: they 
model evolution through time of (say) a heat distribution. A heat distribution 
is an extensive quantity, and does not necessarily have a density function, 
which is an intensive quantity; the most important of all distributions, the 
point distributions (or Dirac distributions), for instance, do not. For the case 
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of the heat equation, it is well known that the evolution through time of any 
distribution "instantaneously" (i.e., after any positive lapse of time, t > 0) 
leads to distributions that do have smooth density functions. But in SDG, 
we are interested also in what happens after a nilpotent lapse of time. In 
more computational terms, we are interested in the Taylor expansion of the 
solutions of evolution equations. For this, it is necessary to stay within one 
vector space, that of distributions. 

The vector space of "distributions of compact support" on any object M 
can be introduced purely synthetically (see [1^ p. 393, or p. 94) as the 



i?-linear dual of the vector space i?*^ (which internally represents the vector 
space of smooth functions on M). What follows could, to a certain extent 
(in particular for the wave equation), be treated purely synthetically. 

Presently, we shall only be interested in distributions on R, , and R^, 
so for the presentation, we have chosen to assume that we are working in a 
sufficiently good "well-adapted" model E of SDG, containing the category 
of smooth manifolds as a full subcategory. In such models, for any given 
manifold M, we could define the linear subspace T>{M) of i?*^ consisting 
of functions with compact support, (the "test functions"). Then the vector 
space of distributions on M, V{M), is taken to be the i?- linear dual of V{M). 

One could take an alternative, slightly more concrete, approach: namely, 
take a model E of SDG which contains the category of smooth manifolds as 
above, but which also contains the category of Convenient Vector Spaces 0] 
and the smooth maps between them as a full subcategory. The embedding is 
to preserve the cartesian closed structure. Such models do exist: we provided 
in , such an embedding of Convenient Vector Spaces into the "Cahiers" 
topos of Dubuc 0. Note that the usual topological (Frechet) vector spaces of 
smooth functions, test functions, distributions, etc. on a smooth manifold M 
have canonical structure of Convenient Vector Spaces. In such a model, we 
can construct internal functions, say curves f : R ^ V{M), by constructing, 
externally, a function by an "excluded middle" recipe of the form 



f(t) = .. ift^O;/(t) = .. ift = 



and then proving smoothness of / by a usual limit argument. 

We have to resort to this kind of "external" constructions only for the heat 
equation, and there our embedding from 0, is not quite good enough, 
since it does not take manifolds with boundary into account; for the heat 
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equation, one constructs externally an "evolution" map 

i?>o ^ V\R) 

by an excluded middle recipe. 

So, for the justification of our treatment of the heat equation, we need 
an extension (hopefully forthcoming) of our work 0, i.e., we need to 
construct a Cahiers-like topos that includes also manifolds with boundary, 
and then to construct an embedding of Convenient Vector Spaces into that 
"extended" Cahiers Topos. (Maybe even the Cahiers Topos itself will be 
good enough.) 

For what follows about wave equation, the Cahiers Topos, and the em- 
bedding of Convenient Vector Spaces into it, is sufficient; in fact, for these 
equations, a purely synthetic treatment alluded to will be sufficient, since the 
distributions considered there are all of compact support. 

As stressed by Lawvere in [^, distributions should not be thought of 
as generalized functions: functions are intensive quantities, and transform 
contravariantly; distributions are extensive quantities and transform covari- 
antly. For functions, this is the fact that the "space" of functions on M, 
R'^^ is contravariant in M, by elementary cartesian-closed category theory. 
Similarly, the "space" of distributions of compact support on M is a sub- 
space of R^^' (carved out by the R linearity condition), and so for similar 
elementary reasons is covariant in M. We shall write V^{M) for this sub- 
space. The space of functions of compact support on M is only functorial 
with respect to proper smooth maps, (counterimages of compact set required 
to be compact), and so similarly, the space V\M) of all distributions on M 
is covariant functorial only w.r.to proper maps. The formula for covariant 
functorality looks the same for V and P^; let us make it explicit for the T>' 
case. Let / : M ^ iV be a proper map. The map V'{f) : V'{M) V'{N) 
is described by declaring 

<I?'(/)(/i),0>=</i,0o/>, (9) 

where is a distribution on M, and is a test function on A^, (so o / is a 
test function on M, by properness of /). The brackets denote evaluation of 
distributions on test functions. 

We shall also write just /(/i) instead of P'(/)(/i). 

Recall that a distribution /i on M may be mulitplied by any function 
: M — ^ i?, by the recipe 

< 5( ■ /i,</) >=< ■ >, (10) 
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observing that g ■ (p is a. test function (has compact support) if is. 

If X is a vector field on M, one defines the directional derivative DxifJ') 
of a distribution /i on M by the formula 

< Dx{lJ^),(j)>= - < f^,Dx{(j)) > . (11) 

This in particular applies to the vector field d/dx on R, and reads here 
< fi', (p >= — < fi,(j)' > (0' denoting the ordinary derivative of the function 
0). One has the following Leibniz rule: 

Dx{f-f^)=Dxif)-lx + f-Dx{f) (12) 

for any distribution fi and function / on M. This is an elementary con- 
sequence of the Leibniz rule for directional derivatives Dx of functions on 
M. 

Remark. The equation (|TTp becomes a theorem, rather than a definition, if 
one takes the following line of reasoning: let F be a covariant functor from 
microlinear spaces (and invertible maps between them) to Euclidean vector 
spaces. Then one may define the Lie derivative along X, Lx{<y), as a map 
F[M) —>■ F{M). For the functor F = T>', Lx becomes the Dx described. 
We shall not pursue this line further here. 

Applying Dx twice leads to 

< DxiDxifi)), >=< /i, DxiDxm > 
In particular, for fi a distribution on i?" 

< dydxi\ij.),(f)>=< ij,dydxi^{(t)) > 
and therefore for the Laplace operator A = J^d'^/d'^^i = divo grad, we put 

<A(/i),0>=</i,A(0)>. (13) 

The following Proposition is an application of the covariant functorality 
of the functor T>c, which will be used in connection with the wave equation in 
dimension 2. We consider the (orthogonal) projection p : i?^ — > i?^ onto the 
xy-plane. (It is not a proper map, so functorality only works for compactly 
supported distributions.) 
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Proposition 12 For any distribution S (of compact support) on , 

p{A{S)) = A{p{S)). 



(The same result holds for any orthogonal projection p of onto any linear 
subspace; the proof is virtually the same, if one uses invariance of A under 
orthogonal transformations.) 

Proof. Let ip be any test function on R^. Then 

< piA{S)),ij >=< A{S),iIj o p >=< S, A(V^ op)> (14) 

But, with ip = ip{x,y), o p is just ■?/', considered as a function of x,y,z 
which happens not to depend on z; so 

OX oy oz 

the last term vanishes because ip does not depend on z, so the equation 
continues 

OX oy 

So the right hand expression in (|T^) may be rewritten as 

=< S, Aiij) o p >=< p{S), A(^) >=< A{p{S)), ij >, 
from which the result follows. 



3.1 Spheres and balls as distributions 

For a,b G R, we let [a, b] denote the distribution f ^ f{x) dx. Such 
distributions on the line, we of course call intervals; the length of an interval 
[a, b] is defined to be 6 — a. Note that the interval [a, b] as a distribution 
is not quite the same as the order theoretic interval, i.e., the subset of R 
consisting of x with a < x < b. For instance, the order theoretic interval 
from to contains all nilpotent elements, whereas the distribution [0, 0] is 
the zero distribution. The distribution theoretic interval [a, b] contains more 
information about a and b than does the order theoretic one. We consider 
the question to which extent [a, b] determines the endpoints. The answer is 
contained in 
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Proposition 13 Let [ai,6i] and [025^2] be two intervals in the distribution 
theoretic sense. They are equal as distributions if and only if they have same 
length, bi—ai = 62 — 0.2 (= ^ say), andl-{ai — a2) = (this then also implies 
/. (6,-62) = 0). 

Proof. Assume [ai,6i] = [(22,62]. The statement about length follows 
immediately by applying each of these two distributions to the function / 
which is constant 1. Generally, we have for any function / that 



■f)2 

f{x) dx = I f{x) dx 



bi 
ai 



f{t + a2 - fli) dt, 

by making the change of variables t = x + ai — a2. Subtracting, we get 
= / {f{x) - f{x + a2- oi)) dx. 

J ai 

Apply this equation to the function f{x) = x, we get 

= / (x — (x + 02 — ai)) = (ai — 02) ■ (61 — ai) = (ai — 02) ■ /. 



ai 



Conversely, assume 61 — ai = 62 — 02 (= ^5 say), and = / ■ (ai — 02). For any 
function /, we calculate the values of the distribution [ai, 61] on /. We have 



[ai,6i](/) = (61-ai) [ f{a^ + t-{bi-ai))dt = l I f{a^+t-l)dt. 
Similarly 



[a2,62](/) = / C f{a2 + t-l) dt. 



The difference is 



I (\f{a^ + t-l)- f{a2 + t-l))dt. (15) 

JO 

By Hadamard's Lemma, /(ai + t • Z) — /(a2 + t ■ I) may be written as (ai — 
0'2) ■ g{0'i,ci2,t) for some function g, and so the integral (|I5|) can be written 

as 
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= Z • (ai - aa) / g{ai,a2,t) dt, 
Jo 

which vanishes if Z • (oi — 02) = 0. 

The assertions about bi — 62 is similar. 

Note the following Corollaries: First, if the length 61 — ai of an interval 
[ai,6i] is invertible (positive, say), then the endpoints ai, 61 are uniquely 
determined by the distribution [ai, 61]. Secondly, for any ^1,^2, we have 

[-ti,ti] = [-t2,t2] implies ti = ^2- 

In fact, by the Proposition, their lengths must be equal, i.e., 2ti = 2t2- The 
distribution [—t,t] will appear below under the name i?t, "the ball of radius 
t in dimension One" . 

We shall also consider such "balls" in dimension Two and Three, where, 
however, t cannot in general be recovered from the distribution, unless t is 
strictly positive. 

Wc fix a positive integer n. Wc shall consider the sphere St of radius t, 
and the ball Bt of radius t, for any t E R, as distributions on R"^ (of compact 
support, in fact), in the following sense: 




< Bt, -0 >= / i'(x)dx = r / ^jJ{t ■ u) du, 

JBt JBi 

where du refers to the surface element of the unit sphere 5*1 in the first 
equation and to the volume element of the unit ball Bi in the second. The 
expressions involving /^^ and J^^ are to be understood symbolically, unless 
t > 0; ii t > 0, they make sense hterally as integrals over sphere and ball, 
respectively, of radius t, with dx denoting surface-, resp. volume element. 
But the expression on the right in both equations make sense for any t, and 
so the distributions St and Bt are defined for all t; in particular, for nilpotent 
ones. 

It is natural to consider also the following distributions S* and on i?" 
(likewise of compact support): 

< S\ip >= / ip{t ■ u) du. 
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, ip >= / fpit ■ u) du. 



For t > 0, they may, modulo factors of the type 47r, be considered as "average 
over S't" and "average over Bt \ respectively, since S*' differs from St by a 
factor which is just the surface area of St (modulo the factor of type 
47r), and similarly for 5*. 

Note that S^ = Si and — Bi. And also note that the definition of 5** 
and 5* can be formulated as 



where Ht : i?" — > i?" is the homothetic transformation u ^ t ■ u, and where 
we are using the covariant functorality of distributions of compact support. 

For low dimensions, we shall describe the distributions St, Bt, 5** and B^ 
explicitly: 

Dimension 1 



S'^HtiSi) ,B'^Ht{Bi), 



< St,^ >^ 7/j{-t) + 7/j{t) 




<S\iP>^iPi-t)+iPit) 




Dimension 2 





■t r-2-K 

/ %l){s COS 9, ssm. 9) s d9 ds 



< 



S\ijj>^ ^jj {t cos 9, t sin 9) d9 



<B\iIj> 




ip{ts cos 9,t s sin 9) s d9 ds 



Dimension 3 



<St,^l^> 




I / xjj {t cos 9 sin (f),t sin 9 sin (l),t cos (f))t'^ sin (f) d9 d(f) 
Jo 
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< Bt^ -0 >= / / / '0(5 cos 9 sin 0, 5 sin 9 sin 0, 5 cos 4>) sin dcf) ds 
Jo Jo Jo 

<S*,'4>>— j I ■ilj{t cos 9 sin (f),t sin 9 sin (f),t cos (f)) sin (f) d9 d(f) 
Jo Jo 

< B^^i/j >= / / / il>{ts cos 9 sin (p,ts sin 9 sin (p,ts cos (j)) sin (p d9 d(p ds. 
Jo Jo Jo 

Notice that these formulas make sense for all t (positive, negative, nilpo- 
tent, ... ), using the standard convention : = — J^), whereas set-theoreti- 
cally St and Bt (as point sets) only make good sense for t > 0. 

It is clear from the very definition that St — t'^~^S* and Bt — t^B^ (in any 
dimension n); but since we are interested also in t's that are not invertible, 
St and S* cannot be defined in terms of each other. 

Note also that Sq = Bq = 0, whereas S^ and B^ are constants times 
the Dirac distribution at the origin 0. The constants are the "area" of the 
unit sphere, or the "volume" of the unit ball, in the appropriate dimension. 
Explicitly, 

5° = 2 -5(0), 27r-(5(0), 47r-(5(0), 

and 

S° = 2-5(0), 7r.5(0), ^-(5(0) 

in dimensions 1,2, and 3, respectively. 

We shall also have occasion to consider the distribution (of compact sup- 
port) t ■ S** on as well as its projection p(t ■ S^) on the xy-plane (using 
functorality of with respect to the projection map p : R^' R^) . For t > 
(more generally, for t invertible), we can give an explicit integral expression 
for it, but note that since S* and t ■ S^ are defined for all t, then so is p{t ■ 5"*), 
whether or not we have such an integral expression. The integral expression 
(for t > 0) goes under the name of Poisson kernel for the wave equation in 
dimension 2 and may be obtained as follows: using the above expression for 
S^ in dimension 3, we have for a test function if) that only depends on a;, y, 
but not on z that 

<t-S^,'4>>— I I ijj{t COS 9 sin (f),t sin 9 sin (f)) ■ t ■ sin (f) d9 dcf). 
Jo Jo 

We then make the change of variables p — tsincf), (f) — arccosp/i, d(f) — 
dp/ — p^, and then the integral becomes 

"f f'^'^ ■0(p cos 9, p sin 9) p d9 dp 



Jo Jo 
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using the explicit form of the ball distribution Bt in dimension 2, we may 
rewrite the right hand side here as 



so that we have, for t > (or even for t invertible), 

Pit ■ S^) = ■ B,. (16) 

4 Vector Calculus 

The Main Theorem of vector calculus is Stokes' Theorem: /^^ u = du, 
for an (n — l)-form, 7 a suitable n-dimensional figure (with appropriate 
measure on it) and d'-f its geometric boundary. In the synthetic context, the 
theorem holds at least for any singular cubical chain 7 : /" M {I^ the 
n-dimensional coordinate cube), because the theorem may then be reduced 
to the fundamental theorem of calculus, which is the only way integration 
enters in the elementary synthetic context; measure theory not being avail- 



able therein. For an account of Stokes' Theorem in this context, see ||I4 
p. 139. Below, we shall apply the result not only for singular cubes, but also 
for singular boxes, like the usual 7 : [0, 27r] x [0, 1] R^, parametrizing the 
unit disk by polar coordinates, 

'-f{6, r) = {r cos6,r sinO). (17) 

We shall need from vector calculus the Gauss-Ostrogradsky "Divergence The- 
orem" 

F) 



flux of F over d'j = (divergence of 



7 

with F a vector field, for the geometric "figure" 7 = the unit ball in i?". For 
the case of the unit ball in i?", the reduction of the Divergence Theorem to 
Stokes' Theorem is a matter of the differential calculus of vector fields, differ- 
ential forms, inner products etc. (See e.g. [0 p. 204). For the convenience 
of the reader, we recall the case n = 2. 

Given a vector field F{x,y) = {F{x,y), G{x,y)) in , apply Stokes' 
Theorem to the differential form 

uj := —G{x, y)dx + F{x, y)dy 
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for the singular rectangle 7 given by ( ]T7| ) above. Then 

7* (dx) = cos 6dr — r sin 6d6 
7* (dy) = sin 6dr + r cos ^(i6' 
7* (fix Ady) = r {dr A (i6') 

Since (ia; = {dG/dy + dF/dx) dx A dy = div (F) (ix A rfy, then 

7*(c/cj) = div (F) r {dr A dO) 

On the other hand, 

7*cu = (F sin ^ - G cos 0)rfr + (F r cos + G r sin ^) rf^, (18) 
(all F, G, and F to be evaluated ar (r cos r sin ^^)). Therefore 

du = C div (F) r dr d9: 



this is /g^ div (F) dA. On the other hand by Stokes' Theorem du = 
Jq^uj which is a curve integral of the 1-form (|18D around the boundary of 
the rectangle [0, 27r] x [0,1]. This curve integral is a sum of four terms 
corresponding to the four sides of the rectangle. Two of these (corresponding 
to the sides 6 = and 6 = 2n) cancel, and the term corresponding to the 
side where r = vanishes because of the r in r {dr A d6), so only the side 
with r = l,0<^<27r remains, and its contribution is, with the correct 
orientation, 

r^(F(cos^,sin^)cose + G'(cos^,sine)sin^) = / F ■ n ds 
Jo JSi 

where n is the outward unit normal of the unit circle. This expression is 
the flux of F over the unit circle, which thus equals the divergence integral 
calculated above. 

We insert for reference two obvious "change of variables" equations. Re- 
call that Ht : —>■ R^ is the homothetic transformation "multiplying by t" . 
We have, for any vector field F on i?" (viewed, via principal part, as a map 
R"" i?"): 

div {Fo Ht) = t- (div F) o Ht, (19) 

and 

t" / (l)oHt = I <p. (20) 

JBi JBt 

We now combine vector calculus with the calculus of the basic ball- and 
sphere-distributions, as introduced in Section 3, to prove 
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Theorem 14 In (for any n), we have, for any t, 

d 



(JjL 



(A = the Laplace operator). 
Proof. We prove first that 

dt 

In fact, for any test function ip, 

< ■ ijj >= t"-^ ■ 4- / V-M du = / (V^)(tu) ■ u du, 
dt dt Jsi Jsi 

(by differentiating under the integral sign and using the chain rule) 

= ■ flux of ((V^) o Ht) over ^i, 

where Ht : i?" — > i?" is the homothetic transformation "multiplying by t" . 
This, by the Divergence Theorem, may be rewritten as 

/ div {{V^) oHt)=t^ f ( div (V^)) o Ht, 

(using (0)) 

= r ■ / (A^) oHt= [ AiP 

JBi JBt 

(by a standard change of variables, cf. (|^)), so 

=< Bt, Alp >=< r ■ B\ Alp >= r- < AB\ i) > . 

From 

^ d^t^ ^ ^ ^^t^ ^ ^ 

dt 

we may of course conclude the desired equality, by cancelling on both 
sides, if t is invertible; but we want the equation for all t. We can get 
this from "Lavendhomme's principle" , which says that if f : R R satisfies 
t- f{t) = for all t, then f{t) is constantly 0. This principle was derived from 
the integration axiom purely synthetically by Lavendhomme in [|ll|] p. 25. So 
the claim of the Theorem is valid for all t. 
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We collect information about t-derivatives of the four basic distributions 
St , Bf, and -B* in i?". The results are valid for any n and any t. For 
invertible t (say positive t), some of the statements may be simplified by 
multiplying by t~^, but we prefer having formulae which are universally valid. 

Theorem 15 We have in dimension n for all t: 

j^iB.)^St, (21) 

t-j^{St) = {n-l)St + t-A{Bt), (22) 

t. ^(B') = S' -nB\ (23) 

j^{S') = t ■ A{B% (24) 



In dimension 1, we also have 



JtiSt)=ABt (25) 



Proof. Equation (plf ) is an immediate consequence of the fundamental 
theorem of calculus; e.g. for n = 2, consider the explicit formula for Bt 
given above in Section 3 ("Spheres and balls as distributions"). With /q as 
the outer integral, the d/dt of it is just the inner integral, i.e., exactly the 
exhibited formula (idem) for St- 

For (^), we t-differentiate the equation St = t""^ ■ 5** by the Leibniz rule 
and get (n - 1) ■ • 5* + T"^ • d/dt{S*); so by Theorem |n|, 

d/dt{St) = (n - 1) ■ . S' + ■ t ■ A(5*). 
If we multiply this equation by t, we get 

t ■ d/dt{St) = {n-l)- ■ S' + ■ A(5*); 

using St = t"'~^S^ and Bt = t"'B^, the result follows (note that A commutes 
with multiplication by t). 
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The proof of ( P^D is similar: t-differentiating t"^ ■ = Bt, we get 

n ■ t""^ -B^ + r ■ d/dtB^ = d/dtBt = St, 

(using (^)), so using St = f^'^S^, this equation may be rewritten as 

-{n-B' + t- d/dtB') = ■ S' 

The result now follows by cancelletion of the factor by Lavendhomme's 
principle, and rearranging. 

Next, (l2^ ) is identical to Theorem and is included again for complete- 



ness' sake. 

Finally, (^) follows from p^): the first term vanishes, since n — 1 = 0, 
and in the remaining equation, we may cancel the factor t by Lavendhomme's 
principle. Alternatively, (p5|) can be proved directly, by a very simple calcu- 
lation. 



5 Wave equation 

Let A denote the Laplace operator J29'^/9x'^ on i?". We shall consider the 
wave equation (WE) in _R", (for n = 1, 2, 3), 

= (26) 

as a second order ordinary differential equation on the Euclidean vector space 
V'^{R"-) of distributions of compact support; in other words, we are looking 
for functions 

so that for all t G i?, Q{t) = A{Q{t)) (viewing A as a map V^iR"") 
V'^{R^). We shall only be looking for particular solutions, in fact, so called 
fundamental solutions: solutions whose initial value and initial speed is either 
the Dirac distribution at 0, or 0. Given any other initial value and speed 
— these being both assumed to be distributions of compact support — , the 
corresponding particular solution may, as is well known, be obtained from the 
fundamental solution just by convolution * with these fundamental solutions. 
This follows purely formally from the rules for convolution of distributions 
P and Q, such as Q * 5(0) = Q, D{P * Q) = D{P) * Q, where D is any 
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differential operator on i?" with constant coefficients; and from linearity of 
the convolution, implying that d/dt{Pt * Q) = {d/dt Pt) * Q; see e.g. |ll6| , 
Ch. 3. 

Dimension 1 

Theorem 16 The function R T^'dR) given by 

is a solution of the WE in dimension 1; its initial value and speed are, re- 
spectively (5(0) and 0. 

The function R l^'dR) given by 

t ^ l/2Bt 

is a solution of the WE; its initial value and speed are, respectively, and 

Proof. The statements about the initial values are immediate from the 
explicit integral formulas for Bt and 5** (putting t = 0). The statements 
about the initial speeds are equally immediate from the following formulas 
(13) and (ID for the t-derivatives, (putting t = 0). We have by (H) 



j^{S')=t-AiB'), (27) 
and so by further t differentiation 

^(S') = A(B') + *4(A(B.)); 

now, d/dt and A commute, so we may continue 

= A(E*) + A(t ■ ^5*) = A(5*) + A(5* - 1 • B'), 

using ( p3D with n = 1. Now by linearity of A, the terms involving 5* in the 
last expression cancel, and we are left with 

^(5*) = HS'), (28) 
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which estabhshes WE for and hence also for 1/2 ■ S^. 
Also, by (pID, we have that 



^W) = S., (29) 
and so by further t differentiation 

^(S,) = |(S,) = ATO, 

using (pSD , which establishes WE for Bt and hence for 1/2- Bt. So the theorem 
is proved. 

Dimension 3 

Theorem 17 The function R —>■ VJ^B?) given by 

in 

is a solution of the WE in dimension 3; its initial value and speed are, re- 
spectively, and6{0). 

The function R —>■ V^{R^) given by 

t^ — -(S' + t^ -AiB')) 
An 

is a solution of the WE; its initial value and speed are, respectively, S{0) and 
0. 



Proof. We calculate first d/dt of t ■ S"*, using (p^ : 

^{t-S') = S' + t^-A{B'), (30) 

and so by Theorem (= (12^)), 

^(t-S')=t- A(B') + 2-t- A(B') + ■ A(^B') 
dt'^ dt 

= 3-t- A{B')+t- A{t- j^B') 
= 3-t-A{B')+t-A{S'-3B'), 
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using (P^l), and now by linearity of A, the terms involving A (5*) cancel, so 
we are left with the equation 

^(t-^*) = A(t-5*), (31) 

which establishes WE for t ■ S** and hence for l/47r ■ t ■ S^. The statements 
about initial value and speed are immediate (using ( pOD for the speed). 

Because /dt^ and A commute, it is clear that if 1 1— >■ Q{t) is a distribu- 
tional solution of WE, then so is t i— d/dtQ{t). So since t ■ S"* is a solution, 
then so is its t-derivative (calculated in (|30D above), i.e. S"* + ■ A(_B*) is a 
solution. Its initial value and its initial speed can be found by putting t = 
in (|3lD (note 6 commutes with multiplication by t). 

Dimension 2 

Recall that we considered the orthogonal projection p : ^ R^. Apply- 
ing covariant functorality, we get for any distribution Q on R^ of compact 
support a distribution p{Q) on i?^, also of compact support. 

Theorem 18 The function R VJ^R^) given by 

t^ — -p(t- S') 

47r ' 

is a solution of the WE in dimension 2; its initial value and speed are, re- 
spectively, and5{fS). 

The function R -> V'^{R^) given by t ^^ l/An ■ p{S^ + t^ ■ A(S*)) is 
also a solution of the WE in dimension 2; its initial value and speed are, 
respectively, 5(0) and 0. 

Recall that an explicit integral formula for p{t ■ 5**), for t > 0, was given 
above, in ( |T^ ) ("Poisson kernel"). 

Proof. The fact that the two distributions in question are solutions of 
the WE is immediate from the Proposition |T2| ( "p commutes with A" ) and 
from the fact that T>'^{p) '■ T>'^{R^) 'D'^{R^) is linear, and hence commutes 
with formation of d/dt\ also, T>'^{p) sends Dirac distribution at G -R^ to 
Dirac distribution at G -R^, so the initial values and speeds are as claimed. 

The Taylor Series at t = for the solutions given can be calculated 
directly, but they can more easily be obtained from the formal solution given 
in Proposition |^. 
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6 Heat equation. 



In this section we deal with distributions that do not have compact support 
and we only consider the one-dimensional case. We are thus considering 
solutions for the vector field on the Euclidean vector space Ty{R), whose 
principal part is given by A : T>'[R) — > T>'{R). We consider the particular 
solution K : R>q V{R) whose initial value is the distribution 5(0). Thus, 
referring to the general treatment of solutions for (differential equations given 
by) vector fields, we are considering R — R>q; for the heat equation, one 
cannot do better, as is well known. Also, as mentioned above, we rely on 
external (classical) calculus; namely, we consider the classical "heat kernel" 
function, i.e., the function K : R>q V'{R) given by 



1 

e" 4r for t > 



47rt . . 

(32) 

6{0) otherwise . 



Here, for the case t > 0, we described a function rather than a distribu- 
tion, so here we do make the identification of functions g{x) with distributions 
I— > S'^^ g{x)(f){x) dx. Differentiation of distributions reduces to differentia- 
tion of the representing functions. For t > 0, wc thus have Kt{x) — K(t,x), 
a smooth function in two variables, described by the above expression. It 
satisfies the heat equation 

dK _ d'^K 
dt dx^ ' 

for i > 0. Also the following hmit expression is classical: 

/■oo 
K{t,x)(l){x) dx = (l){0) (33) 
-oo 

for any test function 0. More generally. 

Proposition 19 For any integer n >0, and any test function (p 

lim,_o+ g^J_^K{t,x)<P{x) dx = 0(^")(O). (34) 
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Proof. The case n = is just (P^D; the general case follows by iteration. 
Let us do the case n = 1. Then 

d d 

— J K{t,x)(f){x) dx = J —K(t,x)(j){x)dx 

K{t, x)(j){x) dx 



dx"^ 



(by the heat equation for K) 

/•no 

K{t,x)(p'^^\x) dx 



(by integration by parts.) 

We then use (|33[) , for the test function 0*^^) to conclude (|3^) for n = 1. 



Proposition 20 The function K : i?>o V{R) is smooth. 

Here, smoothness is taken in the following sense (appropriate for conve- 
nient vector spaces): for each test function 0, the function R>o —>■ R given 
by t \-^< K{t), (j) > is smooth. 

Proof. It suffices to prove that K is infinitely often differentiable at 0, since 
smoothness for t > is clear. For fixed t > 0, we let Kt denote the function 
in X described in (the first clause in) (^) above. Thus, < K{t), > is given 
by the integral 

/•oo 1 

e 4t . (f)(^x) dx. (35) 



3 V47rt 

We first notice that, by Hadamard's Lemma, (f){x) = 0(0) + xtplx). By 
linearity, < Kt,(j) >=< Kt, (f)(0) > + < Kt,xilj{x) > . But < ^^4,0(0) >= 
0(0) and this implies that the derivative of < Kt, > at is 

hm (l/t) < Kt, xip (x) > (36) 



To compute this limit, we use the formulas and notations in Lang's book |T0 
with the exception that we use J-' for the Fourier transform. We also use the 
following well known formulae, where all the functions under considerations 
belong to the class S of fast decreasing functions and thus works with no 
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limitations. First, for any pair of functions a, f3 in this class, one has the 
"adjointness" formula 



Furthermore 



J^((l/v^)e-*« )(x) = (1/V4^)e~" 
^(x^(a;))(0=^(^W)'(0 

To show the existence of the limit, we compute, using 
and (|8l) 

/oo „ 
-oo 

oo 



(37) 
(38) 
(39) 

), adjointness. 



< Kt, xipi^x) > 



oo 
oo 



oo 
oo 



((l/v^)e-*« mx^{x)md^ 

/oo „ 
(l/v^)e-*« ^Hm)d^, 
-oo 

The last step uses integration by parts. Using (|39D, this may be rewritten as 

/oo 
(l/v^)e-*«V(V^')(0 
-oo 

2t I (l/V4^)e-"'/V(a;) dx 



using adjointness and (|37| ) in the last step. Now we divide by t, as requested 
in (p6|), and let t — > 0+. Using (|33|), we thus get that the limit in ( ^61) equals 



lim 2 



(1/V4^)e-" /V(a;) rfx = 2 ■ z^'(O). 

But since 0(x) = 0(0) + x ■ ipi^), 2 ■ ip'i^) = 0"(O), This proves that the limit 
in (|36D exists and equals 0"(O); we conclude that 



^hm (l/t)[< Kt,<P>- 0(0)] = 0"(O). 



(40) 
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To better understand what has been done and to develop this matter 
further, let us define for every t > 

fit) =<Kt,(P> 
We can summarize the results of this section as follows 

/'(o) = r (0) 

Recall from Proposition ^ that 

=< 7^^,0(2") > 

and thus, by going to the limit when t ^ 0~^, 

hm = 0(2n)(o) 

These results suffice to summarize the result in the present Section in the 
following way: 

Corollary 21 The function f is smooth and, furthermore, /*•"•' (0) = 0*^^")(O). 

Proof. Let us show, for instance, that /"(O) exists and equals 0""(O). 
Using the previous results, - /'(O)] = (l/t)[< i^t,^" > - 0"(O)] 

and this implies the corollary, by going to the limit when t ^ and using 
with the function 6" instead of d). Now, iterate. 



The idea to use Fourier transform to prove smoothness was pointed out 
to us by H. Stetkaer and E. Skibsted. 

Summarizing: we have a smooth function K : i?>o —>■ T>'{R), satisfying 
the heat equation dK/dt{t) = A{K{t)) for all t > 0; for t = 0, this follows 
from Proposition |l^. By the assumed fullness of the embedding of smooth 
manifolds with boundary and convenient vector spaces into the model of 
SDG, we have the desired solution internally in the model. We may then 
ask for the values of K for nilpotent t. The answer can be deduced from 
the Taylor Series at for the function K, and the coefficients can be read 
off from Proposition alternatively, by the uniqueness of formal solutions 
(Theorem 0), they can be read off from the formal solution we know already 
from Proposition |^. In any case, we get for nilpotent t 

Kit) = 6(0) + t ■ A(5(0)) + ^A2(<5(0)) + . . . (41) 
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the series being a finite sum, since t is nilpotent. In particular, for d witfi 
rf2 = 0, we liave K{d) = 5{0) + d ■ A(5(0)), or since A = ( )", 

K{d) = 6{0) + d ■ 6{0)" (42) 

In some sense, the motivation for our study of the heat equation in particular 
was to see how 5(0) evolves in nilpotent lapse t of time and specially for t = o? 
with d^ = 0; the answer is (^) (or more generally (^TD). 

Being an extensive quantity, a distribution like (|42|) should be drawable. 
In fact, it can be exhibited as a finite linear combination of Dirac distributions 
6{a) (= "evaluate at a"). This hinges on: 

Proposition 22 Let h'^ = 0. Then 

■ 6{0)" = 6{-h) - 26{0) + 6{h). 

Proof. It suffices to prove, for an arbitrary test function 0, that h'^ ■ 0"(O) = 
4>{—h) — 20(0) + 4>{h)] now just Taylor expand the two outer terms in the 
sum on the right; the terms of odd degree cancel, the terms of even degree 
(0 and two) give the result. (There is a similar result for higher derivatives 
of 5(0): for = 0, 

n 
i=0 

where {n,i) denotes the binomial coefficient n\/i\{n — This hinges on 
some combinatorics with binomial coefficients, cf. p. 63, Problem 16). 

To make a "drawing" of K{d) where d"^ = 0, we assume that d = for 
some h with h'^ = (we shall not deal here with the question whether this 
can always be done). Then 

K{d) = 6{0) + d-6{0)" = 6{0) + h-h^-6{0)" = 6{0) + h-{{6{-h)-26{0) + 6{h)) 

using (|42|) and (^2]). The drawing one can make of S{x) (as for any discrete 
distribution), is a column diagram: erect a column of heigth 1 at x. The 
distribution above then comes about by removing 2h units from the unit 
column at 0, and placing the small columns of heigth h at —h and h. This is 
the beginning of the diffusion of the Dirac distribution. Several other ways of 
exhibiting K{d) as linear combination of Dirac distributions are also possible. 
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Since X>'(i?") is a microlinear and Euclidean i?-module, and A : 'D'{R^) — 
V[R"') is linear, we may apply the general results of Propositions ^ and ^to 
conclude that the formal solution of the heat equation F{t) = A{F{t)) with 
initial value (the distribution) fi, is the series 

/i + tA(/i) + ty2\A\fi) + t'/3\A^{fi) + .... 

Similarly, the formal solution of the wave equation F{t) = A{F{t)) with 
initial value (the distribution) /i, and initial speed the distribution u is the 
series 

fx + tu + e/2\A{^j) + tV3!A(z/) + tV4!A2(/i) + t"" /b\A\v) .... 

Applying (in the one-variable case, say) these formulas to a test function 
in the variable x and to the distributions 5(0) and 5'(0) we obtain the 
following Maclaurin series for the heat equation 

< F(t), >= 0(0) + t0(")(O) + tV2!0("")(O) + . . . . 

Here () refers to the time derivative, whereas ()' to the space derivative d/dx. 
The variable x has been left unexpressed. There is a similar series for the 
wave equation: 

< F(t), >= 0(0) + t0(O) + tV2!0(")(O) + tV3!0(")(O) + tV4!0""(O) + . . . . 



6.1 Simple Transport 

For the sake of completeness, we also consider the function S : R ^ T)'^{R) 
given by t I— > 6{t), the Dirac distribution at t G i?. This is the "fundamental 



solution" for the equation for "simple transport", cf. e.g. [17 



Proposition 23 The function 6 is the solution for the differential equation 
for "simple transport" , 

with initial value 5(0). 

Proof. For any test function 0, 

^ < 5{t), >= |0(t) = 0'(t) =< 5itr, > . 
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